Electron Spin Resonance of SrCu 2 (B03)2 at High Magnetic Field 
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We calculate the electron spin resonance (ESR) spectra of the quasi-two-dimensional dimer spin 
liquid SrCu2(BOs)2 as a function of magnetic field B. Using the standard Lanczos method, we 
solve a Shastry-Sutherland Hamiltonian with additional Dzyaloshinsky-Moriya (DM) terms which 
are crucial to explain different qualitative aspects of the ESR spectra. In particular, a nearest- 
neighbor DM interaction with a non-zero D z component is required to explain the low frequency 
ESR lines for B\\c. This suggests that crystal symmetry is lowered at low temperatures due to a 
structural phase transition. 

PACS numbers: 75.10.Jm, 75.40.Gb, 75.40.Mg, 75.45. +j, 75.50.Mm 
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I. INTRODUCTION 

The spin liquids are states of matter that occur when 
quantum fluctuations are strong enough to avoid any 
type of magnetic ordering. This leads in general to a non- 
degenerate ground state and a finite gap for the spectrum 
of excitations. SrCu2(BC-3)2 is a quasi-two-dimensional 
spin system with a singlet dimer ground stated. This 
compound is at p [resent the only known realization of the 
Shastry-Sutherland models. In this model the effect of 
the quantum fluctuations is amplified by the geometrical 
frustration of the spin lattice. The low energy excita- 
tions of the ground state are local triplets whose "kinetic 
energy" is small compared to the repulsive triplet-triplet 
interaction. The application of a strong magnetic field in- 
duces a quantum phase transition in which the dimerized 
ground state starts to be populated with triplets. The 
magnetic field plays the role of a chemical potential for 
the triplet quasi-particles. In this scenario it is possible 
to study the competition between the "gaseous" triplet 
phase and the crystallization of triplets ("solid phase"). 
The crystallization of triplets gives rise to magnetization 
plateaus that have been observed in SrCu2(B03)2^- 

The anisotropic spin interactions are in general weak 
but they can have a strong effect on a highly frustrated 
system. In particular, as it was already established in 
Refsi 5 i 6 i 7 i 8 i 9 , the inclusion of the nearest neighbor (nn) 
and next-nearest neighbor (nnn) Dzyaloshinsky-Moriya 
(DM) interactions is required to explain some qualitative 
features of the specific heat and electron spin resonance 
(ESR) experiments in SrCu2 (603)2- However, as it was 
noted recently by Cepas et al££>, a lattice symmetry (re- 
flection in the mirror plane containing the c-axis and one 
dimer followed by a tt rotation arround the dimer bond) 
leads to a zero amplitude for the observed single-triplet 
ESR transitions for B\\c. In addition, a level anti-crossing 
between the ground state and the lowest triplet excitation 
is observed for B ~ 20 T. This level anti-crossing implies 
some mixing between two states with different magne- 
tization M z along the tetragonal c-axis, something that 
cannot be explained within the £7(1) invariant models 



(invariant under rotations around the c-axis) for which 
M z is a good quantum number. 

Recent experiments^ revealed additional quantitative 
and qualitative aspects of the ESR transitions. Besides 
the two non-degenerate one-triplet excitations, various 
types of multiple-triplet bound states forming singlets, 
triplets and quintuplets were identified. These measure- 
ments opened the possibility for a direct comparison be- 
tween the observed family of magnetic excitations and 
the theoretical predictions based on the spin model pro- 
posed for SrCu2(B03)2^- In addition, as it is shown 
in the present paper, they provide indirect information 
about the crystal symmetry at low temperatures and the 
role of the spin-lattice coupling as a function of the ap- 
plied magnetic field B. 

The considerable amplitude of the ESR absorption 
lines that according to the crystal symmetry are not ex- 
pected to be observed poses a challenge for finding an 
adequate explanation. Cepas et al. e proposed a mecha- 
nism based on a dynamically generated DM interaction 
induced by the spin-phonon coupling. However, they did 
not provide any comparison between a calculated ESR 
spectrum based on this mechanism and their experimen- 
tal observation. In this paper we suggest that the crystal 
symmetry is lowered due to a structural phase transi- 
tion that occurs in a low temperature region that has 
not yet been explored with X-rays. As a consequence, a 
non-zero c component of the nearest neighbor DM vec- 
tor appears. In a previous paper we showed that this 
component is required to reproduce the measured spe- 
cific heat at low temperatures and high magnetic fields. 
Here, we show that the same component also explains the 
observed single-triplet ESR lines as well as other qualita- 
tive aspects of the ESR spectra as a function of B\\c and 
B\\a. 



II. 



MODEL HAMILTONIAN 



To describe the present system, we consider the fol- 
lowing Heisenberg Hamiltonian in a magnetic field on a 
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Shastry-Sutherland lattice 2 : 

H s = J^Si-Sj + j'^Si-Sj+^s^B-Si 

<i,j> <U>' i 

+ ]T D-(S i xS j )+ ^ D'.(SixSj). (1) 

Here, and indicate that i and j are nn and 

nnn, respectively. In addition to the Shastry-Sutherland 
Hamiltonian, H s includes DM interactions to nn and 
nnn. The corresponding DM vectors are D and D', re- 
spectively. The arrows indicate that bonds have a par- 
ticular orientation as described in Refi The quanti- 
zation axis z is parallel to the c-axis and x to the a- 
axis. The nnn DM interaction has already been con- 
sidered in previous papers 11 to explain the splitting be- 
tween the two single-triplet excitations observed with 
ESR£*i£, far infra-recU^ and inelastic neutron scattering 
measurements^. The value of the DM interaction ob- 
tained from this splitting is: D' = 2.1 K z. According to 
the crystal symmetry^, only the xy component of D is 
non-zero and perpendicular to the corresponding dimer. 
However, as it is explained below, a non-zero z compo- 
nent of D is required to explain the specific heat and the 
ESR data for finite magnetic fields B. 

For D z — and in two dimensions the relevant space 
group of H s in Eq. is p4gm, with a point group 4mm 
at q = point in the Brillouin zone. However, the a-b 
plane in SrCu2 (603)2 containing dimers is slightly buck- 
led and the (three dimensional) space group P42;m is 
more appropriate. Namely, the associated point group 
42m includes the roto-inversion symmetry ICf (rotation 
through ±90° arround the c-axis followed by inversion 
I) which allows for a different orientation of the DM 
interaction with respect to the point group 4mm, and 
thus results in lowering the ground state energy. In zero 
magnetic field H s is moreover time-reversal invariant so 
that the point group at q = should be enlarged to 
Gm = 42m x {E, 0}, where is the time- reversal oper- 
ator and E is the identity operator. 

III. ENERGY SPECTRUM 

Numerical calculations were done using the standard 
Lanczos technique at zero temperature (T = 0) on a 
tilted square cluster of 20 sites. We start by analyzing the 
full energy spectrum for q = 0. Figs, la and lb show the 
calculated energy spectra as a function of B\\c and B\\a, 
respectively. The vertical axis is in units of frequency 
to facilitate the comparison with the ESR experimental 
results by Nojiri et alm&. With the exception of D' z , the 
parameters of the model, J = 74 K, J' = 0.62 J, D = 
(2.2 K, 2.2 K, 3.7 K), and D' = (0, 0, 2.2 K), are the same 
as the ones used to fit the specific heat data for different 
values of the applied magnetic fieldt 

In both cases, there is a clear distinction between states 
that belong to the continuum in the thermodynamic limit 



and those that will be called "localized" states. The bot- 
tom edge of the continuum appears around 1100 GHz 
above the ground state for B = 0. With increasing B 
this edge drops to ~ 450 GHz at B ~ 20 T and then sat- 
urates as a function of B. Since the DM terms are small 
compared to J and J', the localized states can be classi- 
fied according to their approximate total spin quantum 
numbers, S and S z , in the regime B < 20 T. For B < 6 T 
there are at least four singlet and two triplet states split 
due to finite B and D' z . 
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Figure 1: Energy spectrum relative to the ground state energy 
Eq of the Hamiltonian H 3 (Eq. Q) calculated on a 20-site 
cluster with periodic boundary conditions for a) B\\c and b) 
B\\a. The parameters of the model are: J = 74 K, J' = 0.62 J, 
D = (2.2 K, 2.2 K, 3.7 K) and D' = (0, 0, 2.2 K). 

As it is expected from the Zeeman interaction, for B\\c 
(Fig. [IJi) the energy of the two S z = — 1 triplet states 
(O2 and 0\ in the notation of Ref>i£) decreases linearly 
in the applied field. These two states are separated by an 
energy S ~ 90 GHz and their position agrees well with 
experimental ESR lines observed in RefsAiS. Around 
B c = 20 T, the lower triplet state gets mixed with the 
singlet ground state producing a level anti-crossing that 
is a consequence of the finite values of D x and D y . For 
D z = 0, only the lowest triplet state O2 mixes with the 
singlet ground state. This is because the Oi triplet state 
and the singlet ground state belong to different represen- 
tations of the point group Gm- More precisely, a finite 
magnetic field breaks the 0-invariance and reduces the 
symmetry from the full point group Gm to the (magnetic) 
subgroup 42m 1 »£. The ground state and O2 transform 
according to the irreducible representation Ai of 42m, 
while 0\ transforms according to B2. The main differ- 
ence between these representations is in the fact that a 
roto-inversion with respect to the c-axis gives a sign +1 
for the ground state and O2, while it gives —1 for 0\. 
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The hybridization which is induced by a finite value of 
D z = 3.7 K is too small to be observed with the ESR ex- 
periment in the absence of D x and D^terms. These find- 
ings are in agreement with the experimental data&iS. For 
5 1| a the two triplets are nearly degenerate except above 
B = 20 T. Again, in agreement with the experimental 
dataiii, the two triplet states split around B > 20 T 
because of the different hybridization between each of 
them and the ground state. Note that the effect of this 
hybridization becomes significant when the energy dif- 
ference between the triplets and the singlet ground state 
becomes comparable to DM interaction. 

The agreement with the experiment extends even fur- 
ther. The level anti-crossing of triplet and singlet states 
around 646 GHz and 860 GHz for B\\a is also reproduced 
(see Figs. 4b and 7 in RefiiS). Fig.^i also shows a weak 
level anti-crossing of the 0% S z = — 1 triplet state with 
the first singlet excitation located around 520 GHz. This 
effect seems to be too small to be experimentally observ- 
able. However, a much stronger level anti-crossing of the 
2 triplet with a singlet located around and 600 GHz 
is observed at B ~ 2.5 T. Another strong level anti- 
crossing of the 0\ triplet with a singlet bound state is 
observed near B ~ 4 T and 800 GHz. Although the 
values of the magnetic fields of these anti-crossings are 
in good agreement with the experiment, the calculated 
frequencies deviate from the observed values. A compar- 
ison of the calculated ESR spectra for different cluster 
sizes suggests that this deviation is due to finite-size ef- 
fects. For B\\c, the 0\ triplet produces strong level anti- 
crossings for B ~ 4.5 T, v = 600 GHz and B ~ 2.5 T, 
v = 800 GHz. Experimentally, only the later crossing is 
clearly visible in Fig. 4a of Ref»i£. 



IV. SPIN STRUCTURE FACTOR 

To make account of the frequency and the intensity of 
the ESR lines we need to compute the dynamical spin 
structure factor for q = 0, 



1 f 00 
S^w) = -Re / dt, 
Jo 



'(S^S^O)); [i = x,yorz, (2) 



in the direction perpendicular to the applied magnetic 
field. The method that we used to compute S^(uj) is 
described in Refs>i£ii£. Fig. [21 shows the computed ESR 
spectrum as a function of frequency v = lojlix and the 
external magnetic field B along the c- and the a-axis. We 
use this frequency-field type of diagram to directly com- 
pare with the experimental data obtained by Nojiri et 
ali&. The best agreement with the experiment was found 
by normalizing the calculated intensity in a way 

that its integral over all frequencies at a fixed magnetic 
field equals unity. In all figures presenting the calculated 
ESR spectrum such a normalized intensity is visualized 
by the height of the peaks (in arbitrary units). For B\\c 
(see Fig. we obtain a finite spectral weight for the 
two S z = -1 triplet states 2 and 0\. The S ~ 90 GHz 



splitting between the O2 and 0\ states is a consequence 
of the finite D' z = 2.2 K. The level anti-crossing of the 
lowest O2 triplet with the ground state near critical field 
B c ~ 20 T is due to a finite value of intra-dimer DM inter- 
action D x , D y , since these are the only interactions that 
break the rotational symmetry around the z-axis. These 
terms are also responsible for the level anti-crossings at 
B ~ 4.5 T, v = 600 GHz and B ~ 2.5 T, v = 800 GHz. 
Experimentally, only the later one is clearly visible in 
Fig. EJl of RefiiS. The overall effect of finite values of 
D x , D y can be clearly seen comparing Figs. EJa and|21i 
for which D x = D y — 0. For B||c there is no level anti- 
crossing neither with the singlet ground state, nor with 
excited singlets. 

In principle, the D' xl D' terms could also contribute 
to level anti-crossing of 0\ and (or) O2 states with the 
ground state. However, a closer analytical calculation 
shows, that these terms connect the ground state with a 
state, that consists of a product of two triplet states, the 
one with S z = and the other with S z = 1, located on 
the two perpendicular dimers. This state is consequently 
orthogonal to 0\ and O2 states. Therefore, one does not 
expect level anticrossing in the first order in D' x , D' y . 
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Figure 2: Spin structure factor a) S x (uj = 2nv) for S||c and 
b) S z {lo) for B||a. Parameters of the model are: J = 74 K, 
J' = 0.62J, D = (2.2 K,2.2 K,3.7 K) and D' = (0,0,2.2 K). 
Note that the computed ESR spectrum is normalized as 
fduj S' a1 (cj) = 1 for any fixed B (fi = x, z). 

The finite intensities of 0\ and O2 triplets for B\\c are 
a consequence of a non-zero value of D z which requires 
a lower crystal symmetry than the one observed with X- 
raysii In Fig. ^ we show the calculated ESR spectra for 
D z = 0. The 0\ and O2 triplet lines are not observed 
for B\\c while in B\\a case, the lowest triplet excitations 
are clearly visible. Finite values of D x , D y or D' x , D' y 
do not induce these transitions in the lowest order. The 
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reason is that they mix S z = ±1 states with the ground 
state. The non-zero D z term is therefore the only term 
within the given Hamiltonian Q that leads to 0\ and 
02 transitions for B\\c. 

Comparing results of the model Q with "optimally" 
chosen parameters (see Fig. |2l with the experiment in 
RefiiS reveals a good agreement for the line positions, 
and in some cases even matching of level anti-crossings 
with singlet states. The main disagreement with the ex- 
periment is in line intensities. While on the one hand 
ESR measurements of RefiiS show that for B\\c the 0\ 
line is nearly i3-field independent, O2 line shows rather 
strong field dependence: the intensity of the line increases 
with the applied magnetic field. On the other hand, our 
numerical calculations in Fig. [2^i show a nearly constant 
line intensities for both 0\ and O2 lines. 
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Figure 3: The same as in Fig. [5]except D = (0, 0, 3.7 K). 

We focus now on the case B\\a which is shown in 
Fig. |2h>. The 0\ and O2 triplets are nearly degenerate. 
A level anti-crossing with the ground state is well seen 
around B > 20 T. The intensity of the S z = — 1 lines is 
varying non monotonously with B. This non monotonous 
behavior could be due to a weak level anti-crossing with 
localized singlets at u = 600 GHz observed in Fig. ^3. 
Pronounced level anti-crossings are obtained in the up- 
per triplet branches with S z = 1. In contrast to the 
B\\c case, D z is not the only term that leads to transi- 
tions between the ground state and the excited triplet 
states. In Figs. |3Jd and0J) we show S z (u>) for B\\a and 
D = (0,0,3.7 K) and D = (2.2 K,2.2 K,0), respec- 
tively. In both cases we see finite intensities of 0\ and 
O2 triplets, however D = (0,0,3.7 K) leads to a smaller 
intensity than D = (2.2 K, 2.2 K,0). Note that even 
though intensities are presented in arbitrary units, the 
scaling of intensities in all figures is identical to allow 
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Figure 4: The same as in Fig. ^except D = (2.2 K, 2.2 K, 0). 



comparison. 



V. CONCLUSIONS 



In summary, the ESR spectra predicted by model JTJ 
reproduce several aspects of the experimental data ob- 
tained by Nojiri et al. 10 for SrCu2(B03)2- In particular, 
for B\\c, the crystal symmetry breaking interaction D z 
is the only term that leads to finite ESR intensities for 
0\ and O2 triplet states. We have tested other possible 
scenarios that could provide finite ESR intensities for the 
low-lying triplet excitations. One possibility is the intro- 
duction of an anisotropic gyromagnetic g-tensor (Zeeman 
term in Eq. with a different orientation for all 4 spins 
in the unit celUS. In order to get a finite ESR line for B\\c, 
the external field coupled to the g-tensor would have to 
induce a staggered field along the magnetic z-axis, such 
that each spin in a dimer would feel different field ori- 
entation. However, due to the particular structure of 
the g-tensor, which is a consequence of the buckling of 
the ab planes in SrCu 2 (B03)2^, field B\\c only induces 
staggered field component along the x- and y-axis. A 
second possibility is the existence of a small finite angle 
9 between the crystallographic c-axis and the direction 
of the applied magnetic field B due to an error in the 
orientation of the crystal. Taking into account that the 
off-diagonal component of the g-tensor is at most of the 
order of g s ~ 0.05g ~ 0.1, we found almost no detectable 
signal for 8 < 5°. 

The inclusion of the DM interaction D z provides the 
simplest way to explain some qualitative aspects of the 
ESR experiments for SrCu2 (603)2- This simple expla- 
nation has very important experimental consequences. 



5 



The existence of nonzero D z suggests that the system 
should undergo a structural phase transition at low tem- 
peratures that lowers the crystal symmetry. In the new 
phase, the planes containing the c-axis and one dimcr 
are no longer mirror planes. In addition, we also expect 
a strong spin-lattice coupling when the 0\ and O2 triplet 
states get close to the singlet ground state. Such a cou- 
pling could contribute to the stabilization of the different 
plateaus that are observed in the magnetization vs. field 
experiments^. 
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